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! Abstract 

PLh ; The probl em of energy and its localization in general relativity is critically re- 

' examined. The Tolman energy integral for the Eddington spinning rod is analyzed 
in detail and evaluated apart from a single term. It is shown that a higher order 

! iteration is required to find its value. Details of techniques to solve mathematically 

' challenging problems of motion with powerful computing resources are provided. 

■ The next phase of following a system from static to dynamic to final quasi-static 
. state is described. 
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■ criterion, surely the problem of energy in general relativity is a very good one indeed: 
researchers continually return to it in spite of the fact that many might regard it as a 
"solved" problem. An important part of the "problem" aspect derives from the difficult 
issue of energy localizability in general relativity. Einstein [23] , Eddington P] and others 
argued that energy is inherently non-localizable in general relativity while others such 
as Bondi j22j dismissed this notion. Through the years, various authors [Uj j22] have 
addressed the problem and recently, a new group of researchers have entered the debate 



1 Introduction 



It has been said that "a good problem shows its worth by fighting back". By that 
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While the locahzabihty issue is of some interest even for static or stationary systems, 
of particular interest and importance is the issue in the case of a non-stationary gravity- 
wave-emitting system. This is due in part to the fact that although gravity waves 
have been researched intensively over many years, there is only sparse observational 
evidence, and that not actually direct, of the very existence of gravity waves. Moreover 
the existence of gravity waves carrying energy would seemingly constitute a necessary 
condition for the viability of a quantum theory of gravity. Efforts to detect gravitational 
waves directly have emerged as a major research activity in recent years. 

Various aspects of gravity waves led the first author 2J to the hypothesis that energy 
in general relativity is localized in the regions of the non-vanishing energy-momentum 
tensor Tj'^ (henceforth the "localization hypothesis"). For example, plane fronted grav- 
itational waves with parallel rays have been shown to be of the Kerr-Schild class and 
for the latter, it has been shown that all components of the pseudotensor vanish glob- 
ally [15 . This is in stark contrast to the truly tensorial corresponding case of plane 
electromagnetic waves that demonstrate a physically unambiguous flux of energy, both 
theoretically and observationally. Since there is a link between the Bondi news func- 
tion and the pseudotensor [24], the same limitations would be drawn from this route 
of analysis. Another example is that of Bonnor's [30] matching the aspherical Szekeres 
dust collapse metrics to a Schwarzschild exterior indicating a lack of energy flux. In 
this regard, it is also interesting to recall remarks by Pirani [4 : "Suppose, for example 
that a Schwarzschild particle is disturbed from static spherical symmetry by an internal 
agency, radiates for some time, and finally is restored to static spherical symmetry. Is 
its total mass necessarily the same as before?" Very recent support for the localization 



2 



hypothesis comes from the work of Bringley who has shown that Bonnor's directed 
beam of radiation is a Kerr-Schild metric. 

Indications contrary to the locahzation hypothesis come from two directions: firstly, 
there is the gravitational geon construct, a bundle of gravitational waves that self- 
gravitate into a spherical ball with an exterior Schwarzschild metric indicating mass 
yet with a spacetime devoid of any energy-momentum tensor or singularities. However, 
we demonstrated that the time scale for the evolution of the geon was necessarily of the 
same order as the period of the constituent waves themselves, in contradiction with the 
assumptions entering into the construct of the geon itself; a perturbation analysis to 
test the requirement of quasi-stability results in a contradiction [7j. Secondly, there are 
the imploding/exploding Brill waves, again claimed to present mass without an energy- 
momentum tensor. However, this construct proceeds from an instant of time symmetry 
and it is at least questionable whether such an assumption is compatible with nonlin- 
early interacting waves. This second source of possible contradiction to the localization 
hypothesis will have to be considered in great detail. In spite of all the years that have 
since passed and all the papers since written, the question remains. As well, it is to be 
noted that the localization of energy in non-gravitational physics is in the region of non- 
vanishing T°°. Hence the localization hypothesis may be viewed as a generalization of 
established physical phenomena and from this vantage point, not justifiably regarded as 
controversial. 

If the localization hypothesis should prove to be correct, it would have fundamental 
consequences. First, it would imply that gravity waves in vacuum (assuming that they 
exist and there are ample reasons to believe that they do) would not be carriers of 
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energy, in conformity with the Kerr-Schild aspect. This notion challenges the very 
meaning that we give to the word "wave" , as a disturbance that carries energy. Second, 
without an energy aspect to the waves, it would be difficult to argue in favour of a 
quantum theory of gravity. In our earlier papers, we noted that while all particles and 
fields exist within spacetime, gravity in essence is spacetime itself, i.e. it is intrinsically 
different conceptually. Thus, it would not be surprising if the role of energy in general 
relativity should prove to be fundamentally different from that in other areas of physics. 

Tolman jTHj had found an expression for the total energy of a stationary system 
including the contribution from gravity as an integral over the region of the energy- 
momentum tensor. We had focused on this in earlier work 12^ and more recently jH] 
considered the time-rate of change of the Tolman integral at one instant as a vehicle to 
test the necessity condition for an energy loss via gravitational waves from the classic 
Eddington spinning rod The idea was that in the course of the evolution from 
initial stationarity to final (at least asymptotic) stationarity at both points at which 
the Tolman integral registers the energy, the Tolman integral should have undergone 
a change if energy were really being emitted. While the application of the time rate 
of change at one particular instant does not test with rigour the necessity condition 
for an energy loss for reasons that will be discussed later, it could possibly serve as a 
useful at least partial indicator as to how the gravitational contribution to the energy 
is being altered. At that time, the indications were that the Tolman integral was not 
changing. In this paper, we first derive a simpler expression for the time-rate of change 
of the Tolman integral. We then re-examine the Tolman integral more critically and 
find that the complete proper analysis is far more complex than originally envisaged. 
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More information regarding the fields than that provided by Eddington and that we had 
used previously is required and this is furnished. In the course of the analysis, we focus 
upon the techniques of integration over the coordinates of the inertial frame and over 
the co-moving frame, using Dirac delta (a time-honoured helpful calculational tool) and 
step functions as required to properly delineate the source region and motion. These 
techniques will hopefully be found to be useful for researchers in the future. 

We find that all of the terms of the form field/Tj^ product cancel and we are left with 
three time derivatives of a trace quadrupole moment to evaluate, much as in traditional 
flux calculations which employ the untraced quadrupole expression [11^ first derived by 
Einstein |22]- However, for the latter, the quadrupole terms are squared and hence are 
only required to low order of accuracy. By contrast, in the present work, the quadrupole 
trace is not squared for the Tolman integral computation and hence this quantity is 
required to higher order of accuracy, to a second iteration of the field equations. While 
this could possibly be examined in future research, it is valuable to assess the broader 
picture. As Bondi, Bonnor, Feynman and others through the years had noted, the 
truly reliable approach would be to examine the energy of a system at initial and final 
stationary states of the system that has an intermediate wave-emitting phase. If there 
is a change, then the existence and extent of an energy loss can be evaluated without 
ambiguity. The technique described above only examines the rate of change of the 
Tolman integral at an instant of the dynamic phase. Actually, this rate should be 
integrated over the lifetime of the dynamic phase for a complete picture. There is an 
added motivation to do so when one considers that the traditional sources studied have 
been assumed to be periodic as in the present problem, and as Papapetrou showed 
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many years ago, periodic sources are incompatible with the Einstein field equations. 
Thus, the assumption of periodicity is at best an approximation which one hopes is 
adequate when one undertakes analysis with the periodic assumption. Indeed if the 
system has developed over a long period, then questions arise regarding the possibility 
that non-hnearities could develop to lead to a different field than is being envisaged. 
Moreover, the start-up and wind-down phases of the motion are not being considered 
and they could have an important contribution. To this end, we are presently in the 
process of analyzing a system with an explicit relatively short evolution in time from a 
static to a dynamic to a quasi-static final state. This has the advantage of providing the 
desirable total history as described and discussed above as well as mitigating any onset 
of non-linear effects due to the relative brevity of the evolution. 

For completeness, we review some of the fundamentals of the energy problem in 
general relativity and the work of Eddington and Tolman in Sec. |21 In Sec. El the 
manner in which we set out to evaluate the time rate of change of the Tolman integral at 
a particular instant is developed and the framework for its application to the Eddington 
spinning rod is discussed. We illustrate the techniques that are employed to take time 
derivatives, comparing the manner of calculating in both co-moving and inertial frames. 
In Sec. 4, the details of the coordinate dependence for the Eddington problem are given. 
We supply the technical details that form the transition from the required time rate of 
change of the formal Tolman integral to the Maple program that is used to compute it in 
Sec. 5, and we determine the time rate of change of the Tolman integral up to one final 
term. In Sec. 6, we provide a summary and concluding discussions regarding the weak 
points of the assumptions leading to the previous results as well as a discriminating test 
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and suggested future directions. 

2 The Energy Problem 

We recall that the generally covariant energy-momentum conservation laws required for 
general relativity JH] 

T^, = 0. (1) 
can be cast into the form of an ordinary vanishing divergence 

^ {V^iTt + t',))=0 (2) 



by the introduction of the energy- momentum pseudotensor [22] • The pseudotensor 
carries the specifically gravitational contribution to energy and momentum. It is a 
complicated expression involving the metric tensor and its first partial derivatives. Thus, 
global conservation laws can be extended into general relativity but at the expense of 
the inclusion of this non-tensorial object. As a result, the determination of a localized 
expression for energy- momentum becomes problematic as the pseudotensor, unlike Tj', 
can be made to vanish at any pre-selected point by the proper choice of coordinates. 
Einstein considered a weak perturbation of Minkowski space 

gik = Vik + hik (3) 
in conjunction with the field equations 

= ^-^Tt. (4) 

With the field perturbation hik, he computed the pseudotensorial Poynting vector flux to 
lowest non- vanishing order in v/c from a bounded system emitting gravitational waves 



using the right hand side of © below as 

^ I (v^(r° + dV = -cj {^t^) dS^ (5) 

where the integral on the left hand side was regarded as the energy E, extended through- 
out space via tg. 

An alternative approach is to proceed with a multipolar expansion of the field as in 
electromagnetism to express the rate of energy loss to lowest order as 

E = -^JQd.J^)d^^ (6) 



45c5 V^^^^ 7 \dt\ 
where Da/3 is the mass quadrupole tensor [llj. As applied to the Eddington spinning 



rod, both methods yield the energy loss rate 

dE 32GPuj^ 



(7) 



dt 5c5 

where I is the moment of inertia, and u is the angular velocity of the rigid (to lowest 
order) Eddington rod. 

Later, Eddington used a local approach analogous to that of radiation damping 
calculations in electromagnetism. With the conservation laws (0) expressed as {T""^ is 
defined as ^—gT°-^) 



he integrated over a region just beyond the confines of the source and using Gauss' 
theorem, found a useful expression for the time rate of change of the integral of Tq as 

Ur^.v^UT'^l-^h^dv. (9) 
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The beauty of is that it enables one to focus upon the energy- momentum tensor 
part of the energy expression without bringing in the non- locahzed pseudotensorial 
components as in (0). Indeed , after an interesting sequence of argumentation in his 
book P , Eddington concluded that T*"^ was the appropriate density of energy in general 
relativity. Rather than the asymptotic field that is required in the fiux calculation of 
(0), the local field is required in calculations using the RHS of ©• Eddington found 
this in the harmonic gauge as 



-4 



dV. 



(10) 



To bring the elements of the retarded integral to a common time, Eddington performed 
a present time expansion as 



T' 



(1-^) 



n=2 



n\ dt 



(11) 



He considered the rod spinning in the x — y plane and evaluated the right hand side of 
Q with the rod lying along the x axis at t = 0. When the three- volume element dV is 
at X and element dV is at x', the distance between them at time t is 



r = V + a;'^ — 2xx' cos out 



(12) 



which is used in conjunction with pij) and (jH]) and t is set to zero after the 

differentiations. 

Simplifications occur because if T"'' is any non-vanishing component at t = 0, T'ab 
will be an even function of time so odd order derivatives vanish and only the series 
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d_ 
di 



T' 



6 dt^ 



1 (f 



(r^bix^ + x'^ - 2xx'cosujt)^ 



(r^hi^^ + x'^ - 2xx'cosujtf^ 



(13) 



+ 



is required. With the rod spinning in the x — y plane, the required transformation is 



x" = t^t, 



x^ = X — X COS ujt + y sin tut, 



x^ = z — z 



(14) 



X = y = —xsmujt + y cosujt 

where the x'-coordinates are co-moving with the rod and the ^^-coordinates are the 
inertial coordinates. This transformation is of sufficient accuracy for the calculations 
that we will require. 

We will delve into greater detail than did Eddington. As the desired calculation is a 
perturbative one, we begin with the expression of the energy-momentum tensor in the 
co-moving frame to the lowest order, 

T'' = a6iy), 

= a^{x^ - a-)5{y) 



and all other components are zero. We can transform this into the inertial coordinates 
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using (fTi|l to get 

2^00 _ jioo 
= u;xf °° 

(15) 

T^^ = -u^xyf^^ + cosut sin cut 

T22 = u;Vf°o + sin2a;tfi^ 

Expressed in co-moving variables, they are 
2^00 _ jiOO 

= -uj{ycosujt + xsmujt)f°° 
= uj{x cos cut — y sin ujt)T^^ 

(16) 

= -c<;2(ycosc<;t + xsinu;t)2r°° + cos2cjtr^^ 
T"*^^ = —uj'^{x cosujt — y sin cut) (y cos tot + x smujt)T^^ + cos cut sin cut T^"^ 

= uj'^{xcosujt - j/sincut)^f °° + sin^cutf 
Eddington used the information he required from ()16p and again he found the result as 
in dZj). 

Eddington chose to regard the integral of Tq as the total energy of a system. In 
the past, we have referred to this integral as the material or "kinetic" energy but this 
is not a truly adequate description because the root of the determinant of the total 
metric tensor is included rather than that of the spatial part of the metric tensor that 
would provide proper three-volume. Thus, some field contribution is included in the 
Eddington expression. In later years, Tolman (see also piT| for a more elegant 
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derivation) showed that for stationary systems, the total energy is actually 

E = I - T:) dV. (17) 

This is particularly attractive in that it is a wholly localized expression for bounded 
distributions and most importantly, it accounts fully for the contribution to the energy 
from the gravitational field. Moreover, by comparison with (jSl), we can identify the 
pseudotensorial gravitational contribution to the energy as deriving from the integral of 
the trace of the stresses/momentum flux densities. It was this portion of the gravitational 
contribution to the energy of a bounded stationary system that Eddington missed, his 
work pre-dating that of Tolman. We have seen its importance in various ways jl2j . 
While we do not have a localized expression for energy in the non- stationary case, it is 
nevertheless of interest to ascertain whether or not this Tolman integral varies during a 
non-stationary phase. This is because it is the Tolman integral that will be the arbiter of 
energy change when the system eventually returns to stationarity. It must be stressed, 
however, that the information that we can glean at present is incomplete because we 
are incapable of measuring the change in the Tolman integral for the complete history 
stretching from the stationary start to the (asymptotically) stationary end. As well, 
there is the assumption of periodicity with the dynamic field emerging from the source 
at a given instant without regard to prior history. These aspects add further uncertainty. 

There are alternative expressions for the mass of an isolated system. The Bondi 
mass ni(u) is evaluated on the null cone. Since the transition to stationarity is 
an asymptotic one, the evaluation of the mass of system that is left behind after the 
dynamic period would entail the additional condition of allowing t to approach infinity 
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using the Bondi expression, adding to the complexity. 

Another choice is that due to Arnowitt, Deser and Misner (ADM) [20]. Their elegant 
formula expresses the mass of a stationary system as an asymptotic spacelike surface 
integral involving the asymptotic metric tensor components. We j^Ij generalized the 
use of the ADM integral for non-stationary systems by treating it as a Poynting vector 
but we now see this as a pseudotensorial flux with the attendant ambiguities that this 
presents. 

In our view, the Tolman integral is particularly attractive in that it is evaluated over 
the domain of the source and hence it yields the desired final state characteristic when 
the energy-momentum tensor no longer varies in time. This is because the evaluation of 
the Tolman mass at that point of the system that remains behind will be the same from 
that time onwards as the waves continue to infinity. What is being evaluated thereby is 
the end-state active gravitational mass of the material source including the contributions 
from its gravitational field and this is precisely what we seek in gauging whether or not 
there has been an energy loss of the system. One is not burdened by conditions on the 
asymptotic field in this approach. 

3 Formalism to Compute the Variation of the Tol- 
man Integral 

By raising an index, we can readily re-express the conservation law (jH)) as |Sj 
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where 



= \r^'K,,g'' + g%Tt (19) 



As before, in what follows, we will designate as the approximately inertial Cartesian 
coordinates. The Gauss theorem in conjunction with (fT8|l gives 

^ j {T^\^ + r^\^)dV = 2 J T^^dV + J {F^x^ + F^x^)dV (20) 

Similarly, we multiply by two x's and integrate to get 

J T'^^x^x^dV = j {T^^x'' + T'^^x^)dV + J F'^x'^x^dV. (21) 

After taking ^ of (jH)), ehminating ^ / {T^^x^ + T^'^x^)dV by using ^ and setting 
(5 = 7 yields 

/ r^^dV =--^ [ T°VxW- [ F^x^dV--^ [ FVxW 
J 2d{x^YJ J 2dx^ J (22) 

The integral of the spatial trace of the energy momentum tensor is required to complete 
the Tolman integral and its time rate of change is readily found as 

^ ^ T^^dV r^'gkydV = I Tfgp,dV + j r^'go^dV + I T'<^g,,,odV 



dx^ 



(23) 



using (jni) and rjap= diagonal(— 1, —1, — 1) 
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The substitution of ((22) into ^ yields 



dx^ J ^ 2 J dx^ 

1 r 

+ -tttW / F^^'x^'dV (24) 



2 9(x0)2 ^ 
+ j rfKkdV + j T'^'hafidV 

This is a simpler form of the required expression than we had found in earlier work ( [HI ) • 
The second to fifth terms on the RHS of (j^^ are of order hik times Tik and hence of the 
familiar form of the RHS of the Eddington equation ©• However the first term on the 
RHS of must be considered more carefully as there is no metric perturbation in it. 
While the untraced version of this term is only required to low order in the traditional 
flux calculations of the past, this is not the case in the present context because we do 
not square this term to find the time rate of change of the Tolman integral. It is required 
to higher accuracy and this will require a second iteration of the field equations. 

It will also be of interest to evaluate the time rate of change of the "material" angular 
momentum jSH! 

dL-f^ d 



dt 

Using (fT^ and (fnijl and Gauss' theorem, this becomes 

dL'^ 



^ J {x^r'^^ - x^r'^^)dv. (25) 



{F'x^ - F'x')dV (26) 



dt 

for the z component of angular momentum. 

As an illustration of the techniques that we will employ, we compute the time deriva- 
tives of the mass quadrupole moment tensor 



Da,3 = I /i(3a;^a;'^5,^ V - <^a/3X^x^ V) dV (27) 
15 



to the lowest order. 

We can perform the calculation in the co-moving frame (as is commonly done with 
electromagnetic radiation analysis) or in the inertial reference frame. The latter has 
time-dependent parts arising from the density yU = a6{z)S{y){H{a — x) + H{a + a;) — 1) 
where H is the step-function. Equation (fTfj) becomes 

/e re ra 
/ / a5{z)6{y) {H{a -x) + H{a + x)- l){3x''x^ - 6apx''x''6f,^) dxdydz 
e J —e J —a 

a5{z)5{—xsmujt + y cos cut) {Sx'^x'^ — ^apx^x^ dxdydz 



a cos uit 



a cos u)t 



for some e > 0, large enough to enclose the Dirac-delta function. The time-dependence 
of the limits of integration arises by virtue of the confinement of the source as viewed 
in the inertial frame. The rod is confined by two factors, the Dirac-delta function, 
5{y cos ujt~x sin cut), that maintains its position along the co-moving x axis and a double- 
step function, {H{a — {x cos ut + y sin ijjt)) + H {a +{x cos ut + y sin ujt)) — l), that specifies 
its length. The double-step function forces the limits of integration to be time- dependent. 
Had we used fixed limits of integration, x G [—a, a] as opposed to x G [—a cos ut, a cos ut] , 
it would have signified that the rod was expanding as it sweeps through the x-axis 
as shown in Figure H The rule for handling the 5(?/)(i?/-integration, as noted in the 
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Figure 1: The figure on the left shows that the limits of integration is x G [—a, a] whereas 
the figure on the right shows the hmits of integration as x G [—a cos ut, a cos ut] . 

Appendix^, allows us to compute the following: 

2 

Du = -(ja^(3 cos^ tot — 1), 

Di2 = aa^sm{2u!t), 
2 

D22 = -(Ta\2 - Scos^cjt). 

Applying this to (jH)), we get the well-known result of ((Tj), completely independent of time 
and hence without the need to perform time-averaging. 

Had we used fixed limits of integration, we would have found that the resulting 
dE/dt would have been time- dependent with an incorrect value dE/dt = —lOPu^ at 
time t = 0. (Furthermore, because of a cos tut factor in the denominator, the time- 
average would have been infinite as to be expected from an infinitely expanding rod as 

^The rule in section rr.il was aimed mainly at t being in the the neighbourhood of i = 0. However it 
will also apply for —tt/2 < uit < tt/2. Furthermore, without loss of generality, one can argue that the 
calculation is similar for 7r/2 < ut < 37r/2, but with the x's and y's interchanging roles. 
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uot — s> 7r/2.) Thus, it is essential to bear in mind at all times that we are dealing with 
time-varying limits of integration. While it would have appeared simpler to perform the 
calculations in the co- moving frame as is commonly done using 

i*e i*e i*a 



where x{t) = xcosut and y{t) = xsinut, this would not serve our purposes in later 
work. The essential reason is that we will need to take multiple derivatives in the inertial 
spatial directions. Had we insisted on working in co-moving coordinates, we would have 
had to convert back to the inertial coordinates to take the spatial derivatives and then 
convert forward to the co-moving coordinates to perform the integration. As seen from 
the conversion from (fTH|l to this can be rather complicated. Hence, we confine our 
calculations to the inertial coordinates.^ 

4 The Extended Eddington Calculation 

Having established the limits of integration, we now re-examine the Eddington calcula- 
tion (jni) as a second example. Eddington had a one-dimensional problem involving only 
time differentiation and hence was able to perform a simple line-integration. However, 
the computation of the time derivative of the Tolman integral (jl7|) is considerably more 
complex. From the Eddington integral, ()24|) must be subtracted and the latter, which 
contains as in (fT^ . has space as well as time derivatives. Thus, one can no longer 
consider it a one-dimensional problem; one must take spatial derivatives in directions 

parallel to and perpendicular to the rod. Furthermore, the time- derivative of the Dirac- 
^The only exception to this is the integration of the source point where y — 0. 
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delta distribution must be accounted for ^ as well as the double-step functions which 
truncates the rod. Clearly, there is a vast requirement difference between Eddington's 
calculation and that of the Tolman integral derivative. 

The calculation is based on where ?7jj=diag(l,-l,-l,-l) is the flat Minkowski metric 
and ^ 



ij 2 li-i ™ 



hij{x,y,z,t) = -4 / / dx'dy'dz' 

J-eJ^J-a^ [l-Vr)r J ret 

where r = R{x, y; x', t) = ^/ {x — x' cos uty + {y — x' sin uit)'^ + z'^. Eddington used the 
present time expansion (jllj) . This integration is performed in co-moving coordinates. 
The only place that y' comes into play is in the 6{y) and hence it is essentially zero after 
integration. For our purpose, we will not require the expansion beyond the S*'' power. 
The complete expression for hij is 



hij{x,y,z,t) = -4 



— e J —e 



R{Xj z; x', t) 

- ^ (r^ix', r, z', t) - \m,T'-m{x\ r, z\ t\ 

n=2 

It should be emphasized that 

Tij{x\y',z',t) =Tij{x,y,z,t)]^ 

is an inertial entity but it is expressed in co-moving variables. Expressed in this form, 

there is no time-dependence in the Dirac-delta or the double-step functions and hence 

■^When Eddington did his calculation in 1922, the Dirac-delta distribution had not been established. 
^It is possible to use ha = —4 f f° P ^.j 2'?'j^ dx'dy'dz' and integrate in 

the inertial coordinate system but the integration would be extremely difficult; we choose to integrate 
with CO- moving variables. 
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the time-derivatives in the expansion are straight-forward. 

5 Maple Calculations 

At this point, we supply the technical details that form the transition from the required 
time rate of change of the formal Tolman integral (the Eddington term minus (j23)) to 
the Maple program that is used to compute it. 

5.1 The Tolman Integrals 

The original Eddington calculation, 

dt J ^ 2 J dt 

can easily be calculated using the following maple code 

Eddington:= ... intxyz(l/2*TUPPER[i , j] *dif f (hlower [i , j] ,t) ) 

The summations have been omitted for clarity. The full Maple code is available from 
the authors. 

There are a few minor difficulties in calculating ()24|1 . As shown in the Appendix, 
the first time-derivative is commutative with the time-dependent integral. Thus any 
terms having a single time derivative can be readily computed without having to be 
concerned about boundary terms. However, such boundary condition difficulties arising 
from higher derivatives do appear in the first and third terms of plj) . 

We use the definition for F°' in (|T^ to expand the second term of ()24|) to render it 
more manageable by Maple. It is 

^ j F"x"c?y = - j''^^^{^T^h,,,^-T^h^^r^ik)x''dxdydz 

20 



where 77*° = —5*° has been used. It is calculated using the following Maple codes: 

intxyz (dif f ( 1/2*TUPPER [i , j ] *dif f (hlower [i , j ] , xi [alpha] ) , t ) *xi [alpha] ) 
intxyz (dif f (-TUPPER [i , j ] *dif f (HUPPER [k , alpha] , xi [ j ] ) 

*eta[i,k] ,t)*xi [alpha]) 



The last terms are 



di 



^ I T'^'h^, dV 



with the equivalent code of 

intxyz (dif f (TUPPER [alpha , k] *hlower [alpha , k] , t ) ) 

For the third term, we require ()32j) because of the occurrence of the double time- 
derivative. 
1 (9^ r 

2Qt2 J ^ X X 

= -aco^fia) + fi-a)) + T T T ^ {^T'k.^o - ^T'^ hfri,k)x''x" dxdydz 

o Ob o — € o ft 

where the boundary terms are evaluated using the function, 

This is calculated using the regular procedure, intxyz () , and the "boundary" procedure, 
intyzO, as seen in the following code: 

intxyz (dif f ( 1/4*TUPPER [i , j ] *dif f (hlower [i , j ] , t ) * (x~2+y" 2+z~2) 

-1/2*TUPPER [i , j] *dif f (HUPPER [k , 0] , xi [ j ] ) *eta [i , k] * (x"2+y~2+z"2) , t$2) ) 

intyz (1/4*TUPPER [i , j] *dif f (hlower [i , j] , t ) * (x"2+y~2+z"2) 

-1/2*TUPPER [i , j] *dif f (HUPPER [k , 0] , xi [ j ] ) *eta [i , k] * (x"2+y~2+z"2) ) 
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The intyzO lacks x-integration but evaluates the resulting functions at x = —a and 
X = a. It also accounts for the —au'^ factor as seen in the Appendix. 

We find that apart from the quadrupole-like term, the second spatial trace part of 
the Tolman integral has a vanishing time-rate of change after all the remaining terms 
are summed. In the absence of that single term, one would conclude that the Tolman 
integral changes at the same rate as the first part, i.e. as the Eddington value (0), Q- 
A goal in future work would be to evaluate this term by a further iteration of the field 
equations. However, this appears to be a very difficult task and the choice of a different 
dynamical system as we discuss below might obviate the necessity of a second iteration. 

It is to be noted that by breaking up the left hand side of Q into a term that contains 
the determinant of the spatial part of the metric plus a term that carries the remainder 
which is of the form field/Tj^^ product, we find that the time derivative of the remainder 
term is zero. Hence Eddington was actually calculating the time rate of change of what 
should be called the "material" energy, i.e. that which arises from the proper 3-volume 
integral. 

5.2 The Angular-Momentum Integrals 

At this point, we gather additional information from the consideration of angular mo- 
mentum change of the source. Expanded, the expression for dL^'^/dt fl26p is 




The equivalent codes are 
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intxyz (-x/2*TUPPER [i , j ] *dif f (hlower [i , j ] , y) 

-x*dif f (HUPPER [i , 2] , xi [ j ] ) *TUPPER [p , j ] *eta [p , i] ) 
-intxyz (-y/2*TUPPER [i , j ] *dif f (hlower [i , j ] , x) 

-y*dif f (HUPPER [i , 1] , xi [ j ] ) *TUPPER [p , j ] *eta [p , i] ) 

When this is evaluated, we find that 



(28) 



dt 5 

Again, it is worth noting that a spht of the right hand side of the above equation 
as we discussed in the case of the Eddington integral, leads to the analogous result: the 
calculation is the same as for the time rate of change of what we should rightfully call 
the "material" angular momentum. 

Thus we see that ^—^^ as expected for a rigid rod ( ^ = 0) with kinetic energy 
E=^ and angular momentum L=Iui. We had also noted this result in our very early 
work that was based upon pseudotensorial fluxes. However, this relation arises from 
the classical kinetic energy and angular momentum connections with mass, moment 
of inertia and velocities. It is unclear whether or not this lowest order connection is 
appropriate to make the deductions regarding essentially higher order quantities and 
effects. 

In addition, there are important caveats that we consider in the following section. 

6 Summary and Concluding Discussions 

We began by tracing the essential features that have made the issue of energy and its 
localization in general relativity so challenging. The work of Einstein and Eddington 
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were emphasized and as in earlier papers, we focused upon the problem of the Eddington 
spinning rod. An important attractive feature of this source of gravity waves is that to 
lowest order, it can be taken as rigid and hence the time variation derives solely from 
its variation in position. We noted that while Eddington considered the integral of 
as the energy, some years later Tolman showed that for stationary systems, the energy 
is actually the integral of — where the second term accounts for contributions to 
the energy stemming from the presence of the gravitational field. Since the issue under 
consideration is that of the possible gravitational energy loss and since, if it exists, it is 
so incredibly minute in magnitude at least for the kind of sources that we are considering 
here, it is clear that this term is of considerable interest. 

In earlier work, we had noted that various researchers through the years had presented 
the change in the energy as calculated from an initial stationary state to a final (at least 
asymptotically) stationary state as the ultimate criterion for the existence and measure 
of an energy loss. We then argued that a necessary condition for an energy loss was 
therefore a realization of a variation of the Tolman integral during an intervening non- 
stationary period of motion. We had previously used the field information supplied 
by Eddington to compute the change in the Tolman integral at an instant but this 
was inadequate. In this paper, we presented a complete and detailed treatment of the 
fields and found that the combination of all of the Held /Tik product terms produced no 
net time variation of the spatial trace part of the Tolman integral. The final term in 
the equation giving the time rate of change of this trace part of the Tolman formula 
is precisely the trace of the quadrupole moment expression that is used untraced in 
linearized theory with pseudotensor to deduce the "quadrupole formula". However, in 
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the latter treatment, the expression is squared and hence its value is required only to 
lowest order. In the present treatment with the Tolman integral, it is not squared and 
the expression is required to sixth order in v/c. While one can impose the condition of 
rigidity to lowest order as Eddington had done, the field equations impose constraints 
on what is possible for the energy-momentum tensor at higher order and the properties 
of the material come into consideration. We faced this issue many years ago pS] for 
the problem of axially-symmetric free-fall using the pseudotensor and now we see the 
issue of structure appearing anew. The exploration of this issue in the present context 
presents a new and likely difficult challenge that we leave for future study. 
As well, there are other issues to be considered: 

1. Is it sensible to he speaking of a localized energy in the first place? 

Various researchers through the years have denied the possibility of a localized energy 
in general relativity, primarily due to the non-tensorial aspect of the pseudotensor, but 
others such as Bondi |2Z1 argued that energy must be localizable. Moreover, even the 
deniers tacitly admit a certain degree of localization in any case as they speak of a 
material system losing energy and that this loss is due to the gravity waves carrying 
away this energy. In other words, there is an energy to speak of within the material 
source as well as within the waves. 

2. What is the energy of the system during the particularly interesting gravity-wave 
emitting phase? 

The standard derivations of the Tolman integral expression for energy break down 
when the metric is time-dependent. Therefore the energy at that point is unknown. In 
fact it is well to consider what criteria one would employ to decide upon the correctness 
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of any energy expression candidate in the dynamic state. There very well may be an 
additional part beyond the Tolman integral to express the energy for a dynamical system 
that is a localized expression whose density involves the accelerations of the matter. If 
so, it would not always have a non-zero integral: the Tolman integral alone describes the 
total energy for a stationary source such as a rotating disk whose elements, by virtue of 
rotation, do accelerate. A useful goal for future research would be to determine if there 
is an additional part and if so, to deduce its value. 

3. Is the radiation-reaction type of calculation as solidly based in general relativity 
as it is in electromagnetism? 

We regard this as an issue because Maxwell theory is linear and general relativity is 
non-linear. The radiation reaction in Maxwell theory (at least time-averaged) is mean- 
ingfully related to the unambiguous Poynting vector flux of electromagnetic radiation 
energy over even a lengthy period of time, in part because of the linearity. However, 
in the analogous problem in general relativity, there is no assurance that such a calcu- 
lation can yield a meaningful result for anything other than a possibly very short time 
period beyond stationarity during which any non-linearities have not had sufficient time 
to grow to a significant size. Thus, while Eddington jH] argued that ^ m/r for his 
spinning rod was a sufficient condition for a linearized treatment of the energy problem, 
the validity of linearization might be more severely restricted than he had believed. This 
has further ramifications, as we now discuss. 

One might have argued that given the large reservoir of kinetic energy in a rapidly 
rotating Eddington rod as compared to the relatively small amount of gravitational 
contribution to the energy, and given the computed non-zero loss rate of the former (at 
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least within the context of the assumed vahdity of hnearization) , it is inevitable that 
over a long period of time, the loss would accumulate to a large net energy loss value that 
could never be compensated by any gravitational component. However, question 3 above 
is relevant here. It is possible that the kinetic energy loss as computed traditionally holds 
only for a relatively short period beyond stationarity and that during this phase, there 
still remains an equalizing dynamical gravitational component to render the localized 
energy conserved within the source. Moreover, as the non-linearities grow, it is possible 
that the kinetic energy loss component levels out and the vast kinetic energy reservoir 
can never really be tapped. Thus, an experimental test of the localization hypothesis 
would be one of actually observing the equivalent of a rapidly rotating strongly bound 
(i.e. with ^ m/r) system such as an Eddington rod. The rate of change in period 
would have to decrease with the system approaching a state of uniform rotation as in 
classical physics for the viability of the hypothesis. On the other hand, if the decay should 
continue unabated, ultimately draining away the vast reservoir of kinetic energy, then the 
traditional view of an energy loss carried by the gravitational waves would be supported. 
Unfortunately, such strongly bound systems of sufficient strength to yield detectable 
results would be very hard to come by and the far more complex gravitationally bound 
system of a rotating binary such as PSR 1913+16 would not furnish such a test since 
in such a case, the store of gravitational energy is necessarily comparable to that of the 
kinetic energy store, i.e. of the order of m/r. 

We are presently exploring a problem that will hopefully by-pass the difficulties inher- 
ent in the Eddington source. It is based upon an axially symmetric source that evolves 
from a static configuration, evolves into a dynamic phase with an imposed variation in 
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its energy-momentum tensor and then returns after a short time interval to a final rest 
configuration at lowest order. The goal is to integrate the change in the Tolman integral 
over the brief history of evolution to determine whether or not the active gravitational 
mass at the end equals the value at the beginning. By this approach, one avoids the 
uncertainty associated with the potential build-up of non-linearities. Moreover, by com- 
pleting the dynamic phase, the change or constancy in the Tolman integral would answer 
the energy-loss question in a definitive manner. 

7 Appendix 

7.1 Dirac-Delta Integration 

The integral of 



t is in the neighhbourhood oi t — and so for some e > 0, the limits of integrations 
will completely cover the Dirac-delta distribution. Factoring the coefficient out and 
evaluating the integral results in 



where we have assumed that | cos u!t\ = cosut in the neighbourhood of t = 0. Thus 
the rule for integrating 5{—xsmu!t + y cos cut) WRT y is to substitute y — x tan cut and 
divide the entire quantity by cos cut. 




is evaluated using the well-known identity S{ay — h) 



\a\ 



1 



5{y — ^). We assume that 
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7.2 Single Time-Derivative and Integral Commutation 

It is very impractical to integrate before setting time t = 0; setting t = prematurely 
could remove a large number of terms that are required for integration. The identity, 

J -a cos cot J |*=0 J-a"- ) \t=0 

for an arbitrary function f{x,t) is extremely useful. It allows one to bring the time- 
derivative through the time-dependent integral. The proof of this identity is as follows. 
^ Define a function F{x,t) as 

F{x,t) = [ f{x\t)dx' (30) 
Jo 

with the condition that d'^F{x,t) /dxdt = d'^F{x,t)/dtdx. Using this definition in the 
integral on the LHS of (j^^ . we find that 



a cos Lot 



f(x,t)dx = F(x,t)\ -F(x,t)\ 

\x=acosujt \x=— a cos uit 

—acosujt I I 

= F{a COS tut, t) — F{—a cos cut, t). 

When we apply the time- derivative, the chain rule is used on the x = ±a cos urt terms 
and the equation becomes 

d /•'^™«'^* fdFix.t) d , , dFix,t)\ 

OZ J —acosujt ^X OX OZ / \x=acosuit 

/dF(x,t) d , , dF(x,t)\ 

V ox dt Ot / \x=-acosuit (^31) 

= — aoj sin cut (/(a cos cut, t) + /(—a cos cut, t)) 
^,dF{x,t)\ fdF{x,t)- 



dt / \x=acosuit V dt / \x=— a cos LL)t 



^Of course this definition of F{x,t) is not unique. Also, note that we could have defined F{x,t) 
f{x\t)dx' but we cannot be sure that there will be no singularities in the domain [0, x]. 

29 



The link, x = acosut, is applied at the very last step-outside the parenthesis. Hence, 
there is no need to be concerned about time-varying integral limits when we do take the 
time-derivative of (j3U|) . Now, once we set t = 0, we find that the sinut terms go to zero 
simplifying the equation, 

d r=°'^* ^J 1 ( dF{x,t)- 



, Cj i>a COS ujt 

{-/ f{x,t)dx\, = + 



dt / ^x=a,t=0 

dF{x,t)- 



0+ f^, J\ 

Ut ' 3;=— a,t=0 

To show that this is indeed the RHS of ()29|). we refer back to the definition of t) in 
(IHUjl . The expression, 

V dt J \x=a,t=0 \ dt J \x=-a,t=0 \dt J ' 



t=0 



is equivalent to the RHS of (|^ because we can now move the time-derivative through 
the integral. 

7.3 Double Time-Derivative and Integral Commutation 

With a similar approach we can find an expression for the second time-derivative as we 
now demonstrate: 



q2 pa, cos ujt 



idt^ 



f{x, t) dx 

a cos Lijt 



\t=0 

= -auj\fia,0) + fi-a,0))+r\^^^\, dx'. (32) 



dt"^ ) \ t=o 



To prove this, we first define an intermediary function. 



dF{x,t) 

G{x,t) = — — — . (33) 
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When we apply another time-derivative to (j31|) with the help of (|33p. it becomes 

q2 j-a cos ujt Q 

I f{x,t)dx=—i—aujsmujt{f{acosujt,t) + f{—acosLut,t))] 

+ —(G{x,t)\ -G{x,t)\ ). 

(jf \ \x=acosuit \x=—acosLjt/ 

We will not need to expand the first two terms completely because when we set t = 0, 
many parts of 

q2 pa cos Lut . . 

— / /(x, = — ac<j^ coscjt ( /(a coscut, t) + /(— a coscjt, t) j 

'^'^ J -acosujt ^ ^ 

— auj sin i^^^ (^/l*^ cos t) + /(—a cos tut, t) j 
/9G(a;,t)d, , dG{x,t)\ 

(dG{x,t)d, , (9G(x,t) 



a;=a cos 



x=— a cos LDt 



will cancel. With t = 0, the expression simplifies to 

q2 pacosuit 



/acosa;t 
/(x, t) dxl I = - au^ ( /(a, 0) + /(-a, 0) ) 
-acos^^t J |t=0 V J 

dG{x,t)\ ^dG{x,t) 



dt J \ x=a,t=0 \ dt J \x=-a,t=0 

Based on the definition for G{x,t) in (jHHj) . this is 

r pa cos Lut 

'(92F(x,t)\ /d'^F{x,t) 



V (9)!:^ / |x=a,t=0 V dt^ ) U=-a,t=0 

And finally, using (jSO)), we can see that this 

2( , , fd^F{x,t)\ fd^F{x,t)- 

-ac.^(/(a,0) + /(-a,0)) + (^^) |_^,^, " (^^; 

^2 pa 

- I f[x',t}dx' } 

\t=0 



W(/(a,0) + /(-a,0)) + |— / fix',t)dx' 
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is equivalent to the RHS of ((221) • 

7.4 Single Time-Derivative of Step-Function 

Another proof of section 17.21 can be done using step-functions. The basic reason why 
there is time-dependence in the integral limits of (j29|) is because of the step-functions in 
the integrand. An equivalent statement to is 

J {H{a-x) + H{a + x)-l)g{x,y,t)6{y)dydx'j^^^^ 

with the relationship between (j^Hl) and (jH^ being f{x,t) = fl^g{x,y,t)6{y) dy = 
g{x,xtanujt,t)/ cosut. We can move the time derivative through the integral of ()34|1 . 
yielding 

// {H{a — x) + H{a + x) — l)g{x, y, t)5{y) dydx 
-a J —e 

= / {-^yS{x-a)+uyS{x + a))g{x,y,t)S{y) 

J —a J —e 

+ {H{a -x) + H{a + x) - 1)^ (^(2;, y, t)6{y)^ dydx. (35) 

Having now dispensed with the derivative apart from in the very last term, we may set 
t = 0. This implies that x = x, y ^ 0, and H{a ± x) = 1. Thus, 

^1 J {H{a-x) + H{a + x)-l)g{x,y,t)6{y)dydx^^^_^ 



d_ 

dt 



'-ujy6{y)6{x - a) + ujy6{y)6{x + a))g{x, y, 0) 

Using the fact that yS{y) = 0, the first two terms vanish, leading us to (jH^ . 
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7.5 Double Time-Derivative of Step-Function 

The equivalent statement to (jSH) using step functions is 

/ / {H{a — x) + H{a + x) — l)g{x,y,t)5{y) dydx> , 

J -a J-e J |*=0 

-acu2((7(a,0,0) + (7(-a,0,0)) j [^{g{x,y,t)5{y))],dydx. (36) 



To prove this, we first apply another time-derivative to to get 




a J — e 




{H{a — x) + H{a + x) — l)g{x, y, t)6{y) dydx 

2— r/— \ 2— 2 \"^ ^} 2— e/— \ 2 —2 ^{•^ ^) \ / i\r/— \ 

- a) ?/ — - u xd{x + a) -u y —— — - )g{x,y,t)d{y) 

[a — X) —[X + a)/ 

+ 2 (^-uj6{x - a)y + uj5{x + a)y) (5{y)^^^^^^ + y t)) 

+ {H{a -x) + H{a + x) - l)—[g{x, y, t)S{y)j dydx. 

Now, set t = and we get 

/ / {H{a — x) + H{a + x) — l)g{x,y,t)5{y)dydx>. 
J -a J^e J *=o 



Lo'^x5{x — a) — uo^x5{x + a) ]g{x, y, 0)5{y) 




+ 2\^-uj6{x - a)y + uj5{x + a)yj —6{y)g{x, y, 0) 

+ (1) I ^ (^g{x, y, t)5{y)^ } , dyd 



The first four terms of the RHS collapse to —au'^{g{a, 0, 0) + g{—a, 0, 0)) and hence the 
whole expression, 

-auj'^{g{a,0,0) + g{-a,0,0)) + J J ^^^(^g{x,y,t)5{y)^^ ^^_^dydx, 
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is equivalent of the RHS of (jHHj) 
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